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I. INTRODUCTION 


With the current advances of Relativistic Heavy Ion Collider(RHIC) ^ and Cosmology 
in the Laboratory(COSLAB) programs, there is a growing interest in discussing the 
self-organizing nature common to many different physical systems. In particular, the phase 
transition and the defect formation are the paramount examples of the physical phenomena 
due to the self-organizing nature of many-body systems. In the RHIC facilities, vigorous 
analyses are underway to find if the evidences of forming the quark-gluon plasma indeed 
reveal in the experimental data by smashing the heavy ions. Also, the main emphasis 
of the COSLAB program is to exploit the analogies between the phase transitions in the 
universe at the time of the hot Big Bang and the observable transitions in the condensed 
matter systems at low temperatures. Consequently, the physical systems experimented in 
the RHIC and COSLAB programs are very diverse, spanning hadronic matter, superfluid 
Helium, superconductors, liquid crystals, cosmic strings, etc.. Essential theoretical tools to 
study these physical systems are provided by the quantum field theories (QFTs) including 
the quantum chromo^namics (QCD), the toy models of QCD such as Thirring model and 
sine-Gordon model |^, and many other dynamical model theories describing the strongly 
correlated systems. The description of the above self-organizing nature with the QFTs, 
however, cannot rely on the ordinary perturbative approach because the strong correlations 
among the quanta in the formation of new phases are intrinsically nomerturbative going 
beyond the realm of the quantum mechanical von Neumann’s theorem jsf. 

Consideration of phase transitions in early universe has a rather long history Q and 
relevant bubble dynamics has been discussed frequently in the past [^. Recently, the do¬ 
main walls in QCD have been discussed as examples of topological defects which may play 
an inmortant role in the evolution of the early universe soon after the QCD phase transi¬ 
tion Bliiini. In particular, the possibility of primordial magnetic field generation has 
been discussed as a consequence of the existence of QCD domain walls P. Subsequently, a 
claim for the formation of cosmological color magnetic wall in the universe was also made 


by utilizing the QCD instanton solutions |12l. 1131 H 

Due to the topological configurations of the fields, the domain walls and the instanton 
solutions may be stable. Thus, it has been suggested that the QCD domain walls may be de¬ 
tectable at RHIC facilities [^. In the COSLAB, reproducing black hole radiation phenomena 
(he., the so-called Hawking radiation US .) in the table top experiments may be realized by 
the phonon radiation from a dumb hole in acoustic geometry [ig. Also discussed is testing 
various black hole solutions of the Einstein equation in the laboratory Q. Similarly, the 
CO SLAB activities may be applicable as well to study the collisions of the domain walls 
or the so-called bubble collisions in the laboratory. 

Furthermore, the existence of these topological configurations in the universe may give an 
impact on our understanding of the dark matter which has been one of the biggest puzzles in 
the cosmology for more than a decade. Already, the QCD domain wall has been related to 
the axion domain wall . However, the more direct connection between the dark matter and 
the QCD vacuum condensates has been indicated by the recent molar mass estimate of dark 

3 . This 


matter from the strong gravitational lensing in the galaxy cluster 
supported by the analyses of the universal rotation curves in spiral galaxies Q revealed that 
the mass scale of the dark matter particles may be very close to the QCD scale Aqcd 3 i3l- 
Thus, further detailed studies of the QCD vacuum properties such as the QCD domain walls 
and the QCD instanton solutions Bui are necessary to confirm any connection between 
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the QCD vacuum and the dark matter in the universe. 

Motivated by these previous works, we discuss the bubble dynamics in nonabelian Yang- 
Mills theory as a prototype of QCD in this work. In particular, we reanalyze the formation 
of cosmological color magnetic wall with and without following the previous procedure dis¬ 
cussed by Kisslinger, et ah m, in and present both the treacherous points of the 
previous analysis and some new progresses made in the present analysis. Although the 
potential existence of cosmological magnetic wall and the importance of its implication on 
the dark matter and cosmology are undeniable, our analysis indicates that simplihed model 
calculations are not capable of verifying its existence and much more careful studies would 
be necessary to conhrm the possibility. 

The paper is organized as follows. In the next section (Section II), we briefly go over the 
equations of motion and the energy-momentum tensor in the instanton-inspired model of the 
QCD phase transition bubbles. In Section III, we solve the derived equations numerically 
in a plane-wave model describing the collision of the two bubbles in the vicinity of the 
contacting surface: (a) with the Wick rotation procedure followed by the previous analysis 
in Ref. 3, (b) with our own method strictly valid in Lorentzian space. In Section III (a), 
we reproduce the results obtained by the previous analysis and present the error estimates of 
these results to discuss the treacherous points in the analysis. Moreover, we point out that 
the way of solving the system of equations having a constraint should be corrected from the 
previous analysis. By solving the constraint equation hrst and applying it to the dynamical 
one, we show that the numerical results are very different from those in Ref. 141. In Section 
III (b), it is pointed out that the instanton-inspired model studied in Refs. should 

be modified due to the presence of imaginary gluon helds. By constructing a new slightly 
modified Lorentzian model, we present our own results indicating the formation of a bubble 
wall in colliding region. Conclusions and Discussions follow in Section IV. In the appendix, 
more details of error estimates are presented. 


II. INSTANTON-INSPIRED MODEL OF QCD PHASE TRANSITION BUBBLES 


During the time interval between 10“®-10“^s when the universe passed through the critical 
temperature — 150 Mev for the chiral phase transition, the quark-hadron phase transition 
(QHPT) occurs. If QHPT occurs as a hrst order phase transition, bubbles could form. 
Recently, a model for describing bubble collisions has been formulated in Refs. HIQIII. 
Since the main structure of the bubble must be gluonic in nature, pure gluodynamics was 
used. In this section, we briehy review the instanton-inspired model of QCD bubble walls 
in Refs. Bliilli. 

The action for pure glue is 

S = -\j d^xF-F, ( 1 ) 

where 

F^^ = - ig[A^, A^]. ( 2 ) 

For simplicity SU(2) color gauge held can be considered; 


A, = A“(a72) (3) 

with the Pauli matrix cxa satisfying and Tr{a°‘a^) = 25“^. Thus, we have 

F7 = d,Al - d.Al + ge^^^A^^Al, (4) 
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and the field equations in the Lorentz gauge, = 0, become 

d^d^Al + ge’^^\2Ald^Al - A^d.A'^'^) + A^A^^^Al = 0. (5) 

As is well known, the Euclidean equation of Eq. has an instanton solution given by 


= -ft: 


X 


g (x2 + p 2 ) ’ 


( 6 ) 


where p is the instanton size, the 77 “^, symbol is defined in Ref. [21|, and contractions are 
defined in Euclidean metric g^y = diag(l, 1,1,1) for {x'^} = (x*^,x^,x^,x^). Hence, x^ = 
+ x^. Note that the gauge field strength for this solution is self-dual. This instanton 
connects points in two QCD vacua which differ by one unit of topological winding number. 

Inspired by this instanton solution, Johnson, Choi and Kisslinger considered a certain 
class of Euclidean solutions in the form of 

Al{x)=g;^W''{x). (7) 

By substituting it in the Euclidean action of Eq.((T]) and eliminating all fj^^, they get 
1 


1 

'2 


2{d^W,f + {d^W^f + Ag{WWd^W, - W^d^W^) + Zg^W^ 


( 8 ) 

( 9 ) 


in Euclidean space g^i, = diag(l, 1,1,1). Then, by taking the inverse Wick rotation, they 
obtained a Lorentzian Lagrangian corresponding to Eq.© given by 


= - 


2 L 


2{S^W,f + (d,Wp^ + 43(W‘‘W‘'d„W, - 


( 10 ) 


where JE^(x) is the analytic continuation of W^{x) associated with Wick rotation x^ = r = 
it, and contractions are defined in Minkowski space g^^ = diag(—1,1,1,1). Now the field 
equations can be expressed in terms of W^{x) fields as ^ 


d‘^W^ = 2g^W^W^ - 2gW^{d^W^) + 2gW^d^W^, 

and the gauge condition becomes 

= dgW^ - d^Wg. 


(11) 


( 12 ) 


This gauge condition can be expressed as an anti-self-duality condition for the differential 
form of the vector field W, he., *dW = —dW. Further restricting the form of solutions to 


W'^(x) = -x^F(x), 
9 


(13) 


^ It should be pointed out that the direct application of the action principle to the Lagrangian in Eq. cm 
does not reproduce the field equation HU- This follows because, through the ansatz 0, the number of 
unknown functions is reduced from 3 x 4 = 12 for to 4 for so that not all of are independent 
fields. Thus, it should be understood that the instanton-inspired model is defined as the equation dJ. 
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one finds 


(d'^F + AFx ■ dF + 12F^ — 8a; ■ xF^'^ a;^ + 2 (1 — 2a; ■ xF) d^F = 0. (14) 

Here x ■ x = x^ = + x^. The gauge condition of Eq. m becomes 

ei,^apx''d^F = XadpF - X/sd^F. (15) 

Since Eq. dn implies df^F ~ a;^, any solution F of Eq. m satisfies the gauge condition 
above automatically. Multiplying x^ to Eq. m, one obtains 


d‘^F 


-X ■ OF + 12F^ — 8a; • xF^ = 0. 


x^ 


By substituting it into Eq. m again, we have 


d^F- 


X ■ dF \ 

j (1 - 2a; • xF) = 0. 


Thus, one hnds that 


x-dF 

df^F = - Xf,, 

X ■ X 


or F{x) = 


1/2 


2x ■ X —F + X 


^' 


(16) 


(17) 


(18) 


For the hrst case, Eq. (HI becomes equivalent to Eq. (HHl). For the second case, on the other 
hand, Eq. dn becomes 


d^F + AFx ■ dF + 12F^ — 8a; • xF^ = 0, 


(19) 


which is satished by F = (2a; • x)~^. Notice that this equation differs from Eq. (|Tn|l only 
by 4a; ■ dF{F — 1/ (2a; ■ x)) which vanishes in this case. It should be pointed out that, in 
contrast to Ref. [l^, here we have not assumed that F is a function of = a; ■ a; only 
although Eq. da seems to indicate it. One can also see that there is another exact solution 
of Eq. (fTTni given by 

2 - ( 20 ) 

—F + x^ + 

In Ref. □ this solution was considered as an analytic continuation of the Euclidean instan- 
ton solution in Eq.(jni) starting at f = 0. 

From £®*“®(hE) in Eq. (jl()|l . we get the following energy-momentum tensor given by 

a ^ 

= 2d^^W,)d^Wp + g^^dPW^dpW^ - 2g - {W^W, - g^,W^) d^Wp 

-g,uW^W^dpW„^ + 2g^ (gp,W^ - WpW,W^) + 

A \ 

8Fx(^pdi,)F + 2x ■ X dpFduF -|- gpy[12F‘^ -|- AFx ■ dF + x ■ x d^FdpF] 

-8F^ {x ■ X F — 2) {xpXy — X ■ X g^i,) + A {xpXyX ■ dF — x ■ xx{^ydp,)F^ F 


c/2 L 




( 21 ) 
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Note that the energy-momentum tensor used in Ref. Q has wrong sign for the second 
term on the hrst line above. If we insert the exact solution F in Eq. CT into the energy- 
momentum tensor Eq. (EH), appears to be zero. This fact can be expected from the 
self dual property of the instanton solution for the gauge held strength in this system. In 
general, it is well-known that topological gauge helds such as self or anti-self dual solution 
provide the vanishing ener gy d ensity. This is contrary to the results of non-vanishing energy 
density presented in Refs. |l2L IriI|. 

Nevertheless, it should be pointed out that F is singular at points where = 0 

while the corresponding Euclidean solution is regular in whole Euclidean space. Conse¬ 
quently, is not zero everywhere in Minkowski space, but singular at those divergent 
points of F. This singular behavior is analogous to what occurs in type II superconduc¬ 
tivity ji^. Although the magnetic held can not penetrate the superconductor due to the 
Meissner ehect, it is well known that the strong magnetic held can destroy the cooper pairs 
and generate a pinhole in the type II superconductor. In the ideal superconductor this ap¬ 
pears as a singular solution for the magnetic held. However, in realistic superconductor this 
singularity melts due to the interaction with environmental substances. Thus the singular 
energy momentum tensor in Minkowski space might be interpreted to generate a strong 
color magnetic wall in the universe. This sort of interpretation may support the previous 
discussions on the QCD domain walls and the formation of cosmological color magnetic walls 
in Refs. BQ. 

Motivated by this, we would like to discuss the possibility of the formation of this strong 
color magnetic wall in the process of bubble collisions in gluonic dynamics. In the next 
section, we consider two models of bubble collisions and analyze the time evolution of held 
conhguration depending on the initial conditions. 


III. BUBBLE COLLISIONS IN PLANE WAVE MODELS 


In general, the analysis of bubble collision is not an easy task, but involves a lot of 
numerical computations. The authors in Ref. considered a (1 -|- 1) dimensional modeh 
in order to see whether an interior gluonic wall is formed during the collision of two QCD 
bubbles. This model may be appropriate to describe the collision process in the vicinity of 
contact surface, especially, in the case that the size of colliding bubbles is very large. In 
such situations the bubble walls around a colliding surface can be treated as plane waves. 
After we briehy summarize the results of Ref. Q , we present our new analysis. 


A. Instanton-inspired model 

Suppose the collision occurs in the x-direction. Then the gauge helds in the (1-1-1) 
dimensional model Q would be independent of y, ^^-coordinates in the vicinity of contacting 
surface. Taking W2 = W3 = 0 in IT^ in Eq. EH), the authors in Ref. [l^ obtained the 
following equations for Wo{x,t) and Wi{x,t); 

d^Wo = 2g\Wl-W^)Wo-2gWod^W^ + 2gW^dtWi, 


^ This terminology might be misleading since this model is actually defined only in (3-1-1) dimensions due 
to the inappropriateness of 77 -symbols in (1 -|- 1 ) dimensions. 
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d^W, = 2g^(W^ -W^)m + 2glVAW„-2gl'V„d^Wo 


( 22 ) 


with the gauge condition 


d^Wo = dtWi. 

In this model, the authors of Ref. stated that Eq. TO becomes ^ 


d'^F = — d^F - 12F^ + 8(-t^ + x‘^)F^ 

X 


with the gauge constraint 


xdtF + tdxF = 0 . 


With the initial data for the two bubble walls at f = 0 given by 

1 1 


F{x,t)\t=o = 


+ 


{x — 3)2 + p2 (^x + 3)2 + p 2 ’ 


dtF\t=n = 0 


and a periodic boundary condition 


(23) 


(24) 

(25) 


(26) 


F(-10,f) =F(10,f), 


(27) 


the authors in Ref. solved Eq. (El numerically. Since the distance of two peaks is 6, 
we may expect that the time of collision would be f ~ 6 if each peak moves with the speed 
of light. However, due to the smooth prohle of this initial data, the actual collision would 
occur much before this time scale. We have reproduced all of their published results. As 
an example, the function F{x,t) is shown in Fig. ^ which corresponds to their Fig. 2 in 
Ref. [3. Based on this numerical result, they stated that they did hnd a gluonic structure 


evolving at the colliding region although they mentioned the accuracy of the calculation for 
f > 1.0 is limited. 

However, we found that this numerical result is not really reliable for several reasons. 
First of all, notice that the numerical value of F{x,t) drastically changes from the order of 
1 to ~ 10^® around f ~ 1.1. It indicates that some singularities occur around at t ~ 1.1 
so that the accuracy of the numerical calculations is severely limited beyond this point. 
Indeed we hnd that the numerical error estimated by the gauge constraint equation (E3) 
rapidly grows as shown in Fig. El While the gauge condition consists of the two terms, the 
combination of the two terms should not be larger than each individnal term in order for 
the error to be reasonably small. However, Fig. El shows unreasonably large error beyond 
t ~ 1.1. The authors in Ref. have also solved Eq. (1^ numerically for two helds Wq{x^ t) 
and Wi{x,t), and reached the same conclusion of the formation of gluonic walls. Even in 
this case, however, we found severe numerical errors as explained in the appendix. 

Secondly, note that the singular behavior colliding at f ~ 1.1 occurs much before the 
expected collision time scale f ~ 6. As we mentioned earlier this may not be unreasonable 
due to the smooth prohle of the initial data, but in order to clarify whether this indicates 


^ Although we follow their equation to check their numerical results, we found that rewriting Eq. 1221 in 
terms of F{x) gives —instead of in Eq. 12411 . We point out that Eq. is valid only in four 

dimensions. 
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FIG. 1: The evolution of the field F{x,t) in 
the collision of two bubbles. A gluonic wall 
seems to form in the middle. 



FIG. 2: Numerical error estimated by the 
gauge constraint equation, C{x,t) = xdtF + 
tdxF = 0 . 


a violation of causality or not it is necessary to analyze the evolution of the prohle in much 
more detail. 

Finally, we point out that the way of solving the set of coupled equations (1^ and (OK|l as 
above should be corrected. In Ref. F{.x, t) is assumed to be a function oi ^ — only. 

Consequently, the gauge constraint equation (EH) is satished automatically. Now one may 
think that the remaining thing to do is solving the second order partial differential equation 
(1^ for any given initial value of F with vanishing dtF\t=o = 0 as in Eq. (EED- Regardless 
of solving Eq. (El, however, the time evolution of F{x,t) is completely determined by its 
initial value since F{x, t) is assumed to be a function of ^ = x^ — F only. That is, if F = Fq 
at {x,t) = (a;o,0), it should be that, for any position {x,t) satisfying x^ — F = = x^ in 

the future, F{x,t) = F{x^ — F) = F{^q) = Fq. In other words, the initial value of F must 
evolve without any change along the x^ — F = ^ = constant hyperbola. However, we cannot 
observe such behavior in Fig. ^ 

Although F[x, t) is assumed to be a function of ^ only in Ref. 3) show that this as¬ 
sumption is indeed the most general solution of the gauge constraint equation (El. Consider 
the coordinate transformation such that ^ = x'^ — F and r] = tx. Then Eq. (j25ll becomes 


dr,F{Fv) = 0, (28) 

and so the constraint equation can be solved easily, giving that F{x,t) is a function of 
^ = x"^ —F only. Consequently, Eq. El can be rewritten as 

+ 2F'(0 + 3F\0 - = 0, (29) 


where F'{^) = dF{^)/dF This is simply a second order ordinary differential equation with 
respect to F other words, the system of problem which was expressed in the two- 
dimensional space as in Eq. (El is actually in the one-dimensional space described by 
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FIG. 3: The exact solution Eq. (|31 B in FIG. 4: The exact solution Eq. (j3Hl in (x, t)- 

space. space. 


^ : — cx) ~ cx), due to the presence of the gauge constraint equation (EH). Therefore, giving 
the initial data for F{x, f) at f = 0 is equivalent to assigning the value of F{^) for ^ > 0. In 
the case of Eq. (ESI), we have 

^ T7?-qw I 2 + ( 2 

(V^ - 3)^ + (V'C + 3)^ + 

for ^ > 0. Now one can easily check that this function above is not a solution of Eq. ()29j) . 
It implies that one cannot take an arbitrary function of F{x, 0) as an initial data set at 
t = 0 in Eq. ()24j) . The solution space of Eq. alone is much larger than that of Eq. (ESD- 
Only some of them corresponding to a specihc class of initial data satisfy Eq. (EH) as well 
at arbitrary time. 

We have seen that the correct way of solving the coupled equations (El and (El is to 


F{0 




(a) F = —1/2 and F' = —3/8 at ^ = 0. A (b) E = 1 and F' = —3/2 at ^ = 0. Singular 

singular peak occurs at ^ ~ 1.6. peaks occur at ^ ~ —0.8 and —4.0. 

FIG. 5: Some numerical solutions for Eq. (P|) 
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FIG. 6: The time evolution of the numerical solution in Fig. (a) is captured at t = 
0, 1.0, 2.0, and 3.0, respectively. 



FIG. 7: The time evolution of the numerical solution in Fig. [3 (b) is captured at t = 
0, 0.89, 1.5, 1.97, and 2.2, respectively. 


solve Eq. (Ei- What would be the typical behavior of the solutions for Eq. (Ei ? One can 
easily see that 


m 


?,/2 

T 


3/2 


(31) 


is an exact solution of Eq. (1^ . We see that the singular surface propagates along the light 
cone, i.e., ^ = 0 or X = ±t as shown in Fig. EJ When F{^) is not singular at .^ = 0, we may 
have the asymptotic solution near ^ = 0 given by 




(32) 


The numerical solutions for the cases of 6 = 1, —1/2 are shown in Fig. El respectively. In 
(a:, t)-space. Fig. El (a) illustrates a wave initially having two singular peaks at x ~ ±\/T^ = 
±1.26. As shown in Fig. El these two peaks simply propagate away from each other along the 
hyperbolas x^ — F ~ 1.6, and no formation of a wall can be seen in its future evolution since 
\F{x, t) \ monotonically decreases as —oo. For the case of Fig. El (b), the initial wave has 
a single regular peak located at x = 0. This peak quickly grows till t ~ a/TS = 0.894, splits 
into two singular peaks which propagate away along the hyperbolas x^ — ^ ~ —0.8. 

At around t ~ 2 a singular peak is formed again, and splits into two singular peaks which 
propagate away from each other along the hyperbolas x^ — ^ ~ _4 shown in Fig. [7| 

explicitly. 

Therefore, we do not observe the formation of a wall through the collision of two bubbles 
in this model. Even if a wall is formed ‘spontaneously’, as in Fig. (3 it splits into two parts 
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immediately, propagating away from each other faster than the speed of light. The peak does 
not last in the region of collision. This property and the violation of the causality should be 
expected in the region of ^ < 0 basically because F{x,t) is the function of In 

the next section, we discuss where the problematic features of this instanton-inspired model 
come from, and present our new model based on the conventional Lorentzian theory. 


B. New Lorentzian model 


In this section we point out that the instanton-inspired model in which the dynamics of 
bubble collisions were studied above is actually problematic as a Lorentzian theory. Let us 
go back to the construction of such model in Eq. (HH). Under the Wick rotation such as 
x^ = T ^ T = it and a:* —>■ a:* for i = 1, 2, 3 the relationships between helds in Lorentzian 
and Euclidean spaces become 

= r]% = = ■■■ =-i, r]l^ = = 1, (33) 


where is the //-symbol in Lorentzian space corresponding to the Euclidean t’Hooft symbol 
113] • Thus, the relationships between colored vector helds ^4“ in the original Lorentzian 
action in Eq. m and helds in the instanton-inspired model in Eq. m are indeed 


A 


a 


( -iWi -iWt IU 3 -W 2 \ 
-iW 2 -IU 3 -iWt lUi 

V -/IU3 IU 2 -lUi -iWt ) 


(34) 


Since the helds are supposed to be real functions in the instanton-inspired model, it 
implies that the dynamics for the helds in this theory is not usual as in ghost-like 
dynamics in a scalar held system where the scalar held is replaced by an imaginary held, 
e.g., —d^(j)d^(j) = d^'ipd^ip with 0 = ///>. Therefore, although the instanton-inspired model in 
Eq. (fTT|l was constructed in order to describe bubble collisions in Minkowski space starting 
from instantons at t = 0, its dynamics is not the correct Lorentzian one due to the appearance 
of the imaginary electromagnetic helds. Furthermore, it gives rise to the indehniteness in 
the energy density for the helds, e.g., Too ~ —(T“)^ -|- (T“)^. Notice also that the singular 
surface of the solution Eq. (PU|) in the instanton-inspired model expands faster than the 
speed of light. Thus, it appears to violate the causality. 

One may try to redehne some of helds to absorb the imaginary numbers. However, 
such prescription turns out to be impossible as can be immediately seen in Eq. (I34|l . Thus, 
we consider the following model that may describe the Lorentzian bubble dynamics correctly 
in the vicinity of contacting surface; 




^'ijj{x,t) 0 ^ \ 

0 0 y:>{x,t) x{x,t) ■ 

\ 0 0 -x{x,t) J 


This form of helds was obtained from Eq. (I34|l by setting W 2 = IU 3 
symmetry along y and directions. 

With this ansatz Eq. becomes 


(35) 


0 with plane wave 


- TXO - + X^) = 0 (36) 

- 2c/(x(p - (px) - 2///>((p2 + ^ 2 ^ = 0 (37) 

+ 2c/(//^0 - w') - g\{2ip‘^ - + x^) = 0 (38) 

2^(^/^X - XT ) - + X^) = 0, (39) 
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and the Lorentz gauge condition becomes 


%l) — ip' = 


(40) 


Here ip = dtip, ip' = dxip, etc. Thus, we have three second-order partial differential equations 
for three unknown functions (^, V’ and x- Two first-order partial differential equations serve 
as constraint equations. The energy density in this model is given by 


Too{x,t) 


\E%x,t)f + 

2 {gx-ip + dtipf + 2 {-giptlj + dtxf + + x^) 


+Agipxdxip + Agip^^ (gx^ - d^x) + 2 ({d^ipf + {d^xf) 


(41) 


In general, it would be very hard to solve the set of equations fldbHdOj) above. Presumably, 
we have to solve the constraint equations dSi and (HOD first, and then apply it to the 
dynamical equations (jdbHdHj) as was done above for the case of F-dynamics in Eqs. and 
However, we do not know how to solve the constraint equations in general as yet. Let 
us consider some special cases first. 

(i) For the case of y = 0, Eq. (jsni) gives either 93 = 0 or -0 = 0. If 93 = 0, Eqs. (isnD and dHEl) 
are satisfied, Eq. iOD gives -0 = 0, and Eq. (EH) becomes equivalent to a static potential 
problem in abelian theory, he., 'ijj"{x) =0. If "0 = 0, on the other hand, Eq. PH) gives 
ip' = 0, Eqs. (irTf) and (IHH]) are satisfied, and Eq. becomes 

ip{t) + 2g^ip%t) = 0. (42) 

The exact solution for this homogeneous field ip{t) is given by the Jacobi cn function which 
is periodic in time. 



(a) The evolution of the x(x, t) field 


(b) Energy density Too{x,t) 


FIG. 8 : The collision of bubble walls in the simplified case oi p = i/j = 0 and x 7 ^ 0 
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(a) For the case of V’ = 0, Eq. (gni) gives if' = 0 and one finds either x = 0 or (^ = 0 from 
Eq.(jnni)- If X = 0, we have Eq. above. If (^ = 0, on the other hand, Eq. gives 

d‘^x{x,t) - g'^X^{x,t) =0. (43) 


Since the constraint equations are solved already with ip = ip = 0, the remaining thing is to 
solve this single second order partial differential equation only. The numerical evolution of 
two initial peaks for the x(a;, f) held is shown in Fig. |H1 Here the initial data and boundary 
condition are given by 


x{x,t)\t=o = - 


[X 


1 ^ 1 
5)2 + 1 ^ {x + 5)2 + 1 


dtX\t=o = 0, x(-20,t) = x(20,t) (44) 


with g = 1. Although a rather small peak is likely forming in the middle, two bubble walls 
pass away in a short time as in the case of linear waves. Thus, our result seems to indicate 
no formation of a wall through bubble collisions in this simple model. 

{in) For the case of + = 0, Eq. o gives either -0 = 0 or + = 0. If "0 = 0, we have Eq. (14311 
above. If + = 0, on the other hand, Eqs. dSZj) and together with Eq. iOl) result in 

iIj"{x) - 2g‘^il){x)x^{x) = 0, x"{x) - g‘^x{x){x^{x) - 'pj‘^{x)) = 0. (45) 


This case is a static system, which is of no interest for studying dynamical process of bubble 
collisions. Thus, we hnd that the non-vanishing held x{Xi t) with + = -0 = 0 is the simplest 
case for dynamical bubble collisions in the new Lorenzian model (EH), and it appears no 
formation of a bubble wall in such simplihed case. 

Now let us consider bubble collisions in the full model. Since we are unable to solve the 
constraint equations in this general case as yet, we cannot but solve dynamical equations 
only with numerical method. When the coupling constant g is non-vanishing, it can be 
eliminated in Eqs. (I3till4()|) by rescaling helds, he., ip —^ p>lg^ V' V'/fi'; and X ^ x!9- Hence 
we set 5 ^ = 1 in the following analysis. To solve the coupled second order partial differential 
equations ()3bll38j) numerically we take initial data as 


Xo{x) = x{x,t)\t=Q = -— 

{x - 

ipo{x) = dtip\t=o = -1000 


1 ^ 1 
-7)2 + 1 ^ (a;+ 7)2 + 1’ 
tanh(x/100) 
a;2 + 100 


VaM = i’oix) = Xaix), 


(46) 


The initial data for + 0 ( 2 :) and ipo{x) are chosen so that the hrst order equations (ET?|l and 
(gni) are satished. At a; = ±20 we give the periodic boundary conditions 


+(-20, t) = +(20, t), iP{-20, t) = ip{20, t), x(-20, t) = +(20, t). (47) 

Our numerical results for this set of initial data are shown in Fig. El The evolution of the 
gluon helds +, ip and + can be seen in Fig. El (a) in the range of time: f = 0 ~ 8.51, and the 
change of the corresponding energy density is shown in (d). In Fig. El (d) a wall seems to 
appear at around t ~ 4.5. Since the constraint equations (inni) and (gOl) are almost fulhlled 
at least until around t ~ 4.5 as can be checked in Fig. El (b) and (c), this may indicate a 
formation of bubble wall indeed. 

We also have considered some other sets of initial data. Fig. ^1 shows our numerical 
results for +0 = V'o = 0.3+0; 0.05+0; and O.Ol+o with +0 = —100tanh(a;/100)/(a;2 + 100). 
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(a) The evolution of the gluon fields 'tp{x,t), and 





(b) Evaluation of the LHS of 
the constraint equation (EH) 


(c) Evaluation of the LHS of 
the constraint equation dini) 


(d) Energy density Too(x,t) in 
Eq. (1^ with g = 1 


FIG. 9: The collision of bubble walls in the full model 


Formation of bubble wall occurs for all three cases as well although the moment of bubble 
wall formation differs depending on the value of initial data. As in the previous case, the 
hrst-order constraint equations are almost fulhlled until around the time from which bubble 
wall starts to form. Note that, in a longer time, the evolution of two initial peaks for small 
(po and 'ipQ is very much different from that for the case of vanishing ip and 'ijj in Eq. (j4d|l 
even if they are similar in early time. It implies that the presence of all non-vanishing three 
helds </9, -0 and y in our model is necessary for the formation of bubble wall during collisions. 


IV. CONCLUSIONS AND DISCUSSIONS 

In this work, we have considered the collision of gluonic bubbles in the context of the 
instanton-inspired model of QCD phase transition bubbles discussed in Refs. 

In particular, it has been investigated the possibility of bubble wall formation during pure 
gluonic bubble collisions in plane wave approximation. 

Since we are dealing with a system of coupled nonlinear equations of the color dynamics, 
there exist both dynamical equations and constraint equations in this color dynamics. Cor¬ 
rect approach to solve these kinds of system is to solve the constraint equations hrst and 
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(a) ipQ = iPo = 0.3xo 


(b) yjQ = "00 = 0.05xo 


(c) yjQ = "00 = O.Olxo 


FIG. 10: The energy density for various initial values of (p and The formation of bubble wall 
appears at around t ~ 7.5, 10.5, and 11.0, respectively. 


then apply to the dynamical equations. In the case of F-dynamics in Eqs. and (123), 
we can End even an exact solution in this way. Our F-dynamics case corresponds to the 
holonomic constrained system. For other cases, the problems are much harder because the 
constraints are not easily removed. At present, we cannot but solve these problems with 
numerical method. However, we checked that our numerical results satisfy the hrst-order 
constraint equations, up to some small error, till the appearance of bubble wall formation. 

With this correct approach described above we reanalyzed the F-dynamics in Eqs. (12411 
and (123, and found no indication of the bubble wall formation in colliding processes. More¬ 
over, it is pointed out that the instanton-inspired model studied in Refs. |l^ 1^ | (he., 

Eq. with W 2 = H 3 = 0) should be corrected due to the presence of imaginary gluon 
fields, since it leads to the violation of causality. Therefore, we reconsidered the process 
of bubble collisions in a new slightly modihed Lorentzian model (IS3 where we have three 
independent real gluon fields. For the case of y-dynamics with <^9 = -0 = 0 in Eq. (I43|l . 
we did not find any evidence of the bubble wall formation. For more complicated cases in 
Eqs. (in3-(ii3, however, we see some indication of forming the bubble wall. It is likely that 
the presence of all non-vanishing three gluon helds is necessary to have the formation of 
bubble wall in collisions. 

Nevertheless, we stress that much more careful analyses are necessary to confirm our 
results. Although the bubble wall formation occurs in the colliding region of two peaks, 
the location is not exactly the colliding point except for the case of Fig. (a). As can be 
seen in Fig. |21 (b) and (c), the first-order equations are not satished well beyond the time of 
bubble wall appearance. Consequently, the numerical results beyond this time may not be 
trustworthy in Figs. Inland cni The unexpected appearance of peaks at the left end in late 
time evolution in Figs. Inland cni may be related to this increasing error as time goes by. It 
will be of great interest to see how the bubble wall formed actually evolves afterwards. 

Investigations on solving the constraint equations for the complicated cases are in 
progress. The idea of applying canonical transformations to the color-dynamic constrained 
system as in the case of F-dynamics may deserve further consideration. 
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APPENDIX: More Details of Error Estimates 


Since the analysis based on Eq. (USD for bubble collisions is severely constrained both in 
the number of degrees of freedom and the form of the function F{x, t), one may consider an 
extension to the study with the larger degrees of freedom using Wq and Wi. Note that the 
set of equations (E2D and (ESD for Wq and Wi consists of two coupled second-order partial 
differential equations with one first-order equation in time. As explained in the case of 
F-dynamics in Eqs. ()24|1 and dSSD, the solution for Eq. should satisfy the constraint 
equation (OHp as well. In general, not all solutions for Eq. (1^ fulfill Eq. (j2SD- Since it is not 
easy to solve the constraint equation for this case, in contrast to the case of F-dynamics, we 
cannot but solve the dynamical equations ()22|1 numerically for certain initial data. 

Since the evolution equations in (El are second-order, one can take any initial data set 
of Wo{x,0), 1 ^ 1 ( 0 :, 0), dtWo = Wo{x,0) and dtWi = hEi(a;, 0) at f = 0 subject to the first- 
order equation (El evaluated at f = 0 being satisfied (he., lEi(x,0) = O^Wq^x^Q)). The 
case of F-dynamics in Eq. El corresponds to the special case of lEi(a;, 0) = ^a;F(a;,0), 

Ho(x,0) = 1 ^ 1 ( 0 :, 0) = 0, and Wo{x,Q) = —^F(x, 0). In Ref. Q the authors consider a 

somewhat restricted set of initial data such that Wq^x, 0) = c, hEi(x, 0) = w{x), Wo{x, 0) = 0, 
and Wi{x,Q) = 0. This initial data set satisfies the gauge condition Eq. (El, and the 
equations of motion evaluated at f = 0 are written as 

da,w{x) = g{w‘^ - c^) + ^Wo{x,Qi), (48) 

2gc 

dlw{x) = 2g‘^w{nP‘ — c^) -|- lEi(x, 0), (49) 

through which the values of W^i^x^t) and Wi{x,t) at f = 0 can be determined. However, 
the authors in Ref. Q restricted the initial data set further somehow into the case of 
Ho(x, 0) = Wi{x, 0) = 0. For such case the function w{x) is determined as 

w{x) = —ctanh[c 5 f(a: — Xq)]. (50) 


Here Xq is an integration constant. 

As mentioned above, however, one can take in general any function w{x) as an initial 
profile subject to some unknown restriction due to the presence of the unsolved constraint 


equation. Authors in Ref. used a symmetric ansatz for an initial profile of w{x) given 
by 


w{x) = — c {1/2 -|- tanh [cg{x — xq)]} — c (1/2 — tanh [cg{x + xq)]} (51) 


^ The description for the initial data in 



is somewhat misleading as we clarify it in this work. 
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(a) roo(x, 0) (b) -roo(x, t) (c) d^Wo - dtWi 

FIG. 11: A numerical solution up to t = 0.28 in the dynamics of Wq and Wi 


with c = 4, = 1, and a:o = 3. They considered the symmetric bonndary condition, 

ITj(—10,t) = and vanishing initial velocities, dtWi{x,Q) = 0, for i = 0, 1. In this 

plane wave model, the energy-momentnm tensor density TQQ{x,t) is given by 


Too{x,t) = 3{dtWof + - 2dtWod^Wi - {dtWiy - (d^Wo)^ 

+2g \2WldtWo - WoWi{d^Wo + dtWi)] + 2g^Wl{Wl - W^) + (52) 


in terms of the helds Wo{x,t) and Wi{x,t). Here represents the glnonic Lagrangian 
density given by 


+4p [wldtWo + W^d^Wi - WoWi{dtWi + d^Wo)] + 3g\W^ - (53) 

Note that the Lagrangian density in Ref. □ has the opposite sign of it with the absence 
of the hfth term above. 

As can be seen in Fig. [TT] (a), we point ont that the energy prohle corresponding to the 
initial data given above is not in the form of two Inmps that anthors in Ref. 3 expected. 
The nnmerical solntion evolved np to t = 0.28 is shown in Fig. (b). One can see in 
Fig. HI] (c) that the constraint eqnation OxWq — dtWi = 0 is satished well np to t ~ 0.28 
since \dxWQ — dtWi\ is very small in comparison with IclxhFol, l^thFil whose nnmerical valnes 
tnrn ont to be the order of ~ 1. This nnmerical calcnlation was able to be performed np to 
t = 0.28 by using the NDSolve program in MATHEMATICA. However, if we request plotting 
the result further beyond t = 0.28, MATHEMATICA gives it by performing extrapolations 
which are illustrated in Fisf. 1121 

One can see that this extrapolated energy density changes drastically around t = 0.28 
from the order of 10^ to the order of 10^^^, indicating a singular behavior of solution. One 
might think that the^ak shown in Fig.^](a) is an indication of the bubble wall formation 
as discussed in Ref. However, it is hard to think that an initially single peak develops 
a bubble wall. The consideration of other type of initial data giving an energy prohle of 
two peaks would be necessary. Fig. UD (b) also shows that the constraint equation is not 
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(a) -roo(x, t) (b) d^Wo - dtWi (c) Wo{x, t) 

FIG. 12: The extrapolated solution up to t = 5 in the dynamics of IFo and Wi 


satisfied well for this extrapolated solution since becomes the order of 10^® 

which is almost same as the order of \dxWQ — dtWi\. For these reasons, therefore, we think 
that more careful analysis is needed on the bubble collision in the dynamics of Wq and Wi. 
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